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For the classical A'^-vector model, with arbitrary A'', we have computed through order P^'' the high 
temperature expansions of the second field derivative of the susceptibility Xii^^P) on the simple 
cubic and on the body centered cubic lattices. (The A^-vector model is also known as the 0(N) 
symmetric classical spin Heisenberg model or, in quantum field theory, as the lattice 0{N) nonlinear 
sigma model.) By analyzing the expansion of X4(-^i P) on the two lattices, and by carefully allowing 
for the corrections to scaling, we obtain updated estimates of the critical parameters and more 
accurate tests of the hyperscaling relation du{N) + 7(A'^) — 2A4(A'^) = for a range of values of the 
spin dimensionality A^, including N = [the self-avoiding walk model]. A' = 1 [the Ising spin 1/2 
model], N = 2 [the XY model], N = 3 [the classical Heisenberg model]. Using the recently extended 
series for the susceptibility and for the second correlation moment, we also compute the dimensionless 
renormalized four point coupling constants and some universal ratios of scaling correction amplitudes 
in fair agreement with recent renormalization group estimates. 

PACS numbers: 05.50+q, ll.15.Ha, 64.60. Cn, 75.10.Hk 



I. INTRODUCTION 

We have recently extended the computation of high temperature (HT) series for the A^- vector modelEI with arbitrary 
spin dimensionality N on the d-dimensional bipartite lattices, namely on the simple cubic (sc) lattice, on the body 
centered cubic (bcc) lattice and on their d-dimensional generalizations. In previous papers we have tabulated through 
order P'^^ the series for the zero field susceptibility xi^iP) ^-nd for the. second moment of thR-|Correlation function 
H2{N,P) and we have analyzed their critical behavior in the d = 2 casea and in the d = 3 caseEl. Here we present a 
study of the second field derivative of the susceptibility Xii^j P) whose HT expansion on the sc and the bcc lattices we 
have extended through order P^'^ . A study of Xii^, (3) in the d = 2 case had been discussed in Ref.H. It is interesting 
to point out that, in all analyses presented below, the bcc lattice series appear to be better converged than the sc 
lattice series and lead to estimates of critical parameters which are likely to be more accurate. In other words, the 
bcc series seem to be always "effectively longer" □ and therefore give estimates of greater value than the sc series. 

The list of the expansions of X4{N,(3) in d = 3 published up to now is a short one. A decade ago M. Liischer and 
P. Weiszl3(see also Refsfl) derived HT expansions of X4^{N,(3) througii-iiif , for any on the sc lattice in d = 2,3, 
and 4 dimensions by using a linked cluster expansion(LCE) techniqueEnZrEj. In the iV = 1 case, [correspondifiCr-to the 
Ising spin 1/2 model] the series for the. sc. lattice published before our work already extended through /3^Tlj'L3 and 
has been analyzed by various authoraijO. Finally in the Isintt-jaadcl case, a series to order (3^^ on the bcc lattice 
and a series to order on the face centered cubic (fee) latticala'tSliZl have long been available. 

In our calculation we have also used the (vertex renormalized) LCE technique and have developed algorithms 
which are equally efficient in a wide range of space dimensionalities. So far other expansion methods have given 
competitive (or sometimes superior) performance only for discrete site variables and for very simple interactions, on 
two-dimensional or low coordination number lattices. By the LCE method we have produced tables of series expansion 
coefficients given as explicit functions of the spin dimensionality N , with an extension independent of the structure and 
dimensionality of the lattice. Thus we have succeeded in efficiently condensing a large body of information concerning 
infinitely many universality classes. We consider these coefficient tables to be the main result of our work and, in 
spite of their considerable extent, we have reported them in the appendix in order to make each step of our work 
verifiable and reproducible. The size of our computation has been unusually vast: approximately 3x10® topologically 
inequivalent graphs have been listed and evaluated. Nevertheless, we are confident that our series have been correctly 
computed, not only because our codes have been thoroughly tested, but also because N and d enter in the whole 
computational procedure as parameters. As a consequence, at least simple partial checks are available by observing 
that our expansion coefficients, when specialized to iV = 1 agree with tia-series 0(/3"'^^) already available in 3 (as well 
as in 2) dimensions and, for N oo, agree with the spherical modelEEnlj series which can be readily calculated in 
any dimension. More comments on the comparison of our results with the existing series, can be found in our paperS 
devoted to the two-dimensional A^— vector model. 
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A valuable justification of our work is that an increasingly accurate study of the critical behavior of x^i^^P) can 
offer, for all values of N, a sharper test of the hyperscaling exponent relation di'{N) +^{N) — 2A4(7V) = 0. Here j{N) 
and i'{N) characterize the critical singularities in x{N,f3) and ^{N,(3) respectively, while A4(A^) is the gap exponent 
associated with the critical behavior of the higher field derivatives of the free energy. It is also of great interest to 
measure accurately the critical amplitude of X4{N,P), which together with the amplitudes of xi^iP) ^nd (^{N,P), 
enters into the definition of the universal dimensionless renormalized four point coupling constant gr{N). Indeed the 
uncertainties, probably still of the order of 1%, in the value of this quantity might be the main residual source of 
erroiEJO iu the present computation of the critical exponents -aatlun the renormalizatiop-|group (RG) approach by 
the ParisiEa fixed dimension (FD) coupling constant expansiorOcZl. Murray and NickelE2l have recently pushed to 
seven loop order these calculations and the impact of the additional terms on the estimates of the critical exponents 
and of some universal amplitude combinations has been critiC|al|ly assessed by Guida and Zinn-iustinEll. 

As has been stressed many times in the past two decadesESca and, more recently, also in RefB, in order to improve 
the precision of the estimates obtainable from HT expansions not only longer series should be computed, bwjt also 
more careful allowance should be made for the singular corrections to scaling. Their presence is expectedl23 and, 
unsurprisingly, they turn out to be important in various cases. Therefore in this analysis we have also studied their 
role and have estimated their amplitudes in the case of gr{N), both on the sc and the bcc lattice. Moreover, it is of 
some interest to compute the ratios of these correction amplitudes with the analogous quantities for x and ^, which 
define interesting universal quantities, still subject to significant uncertainty and so far not much studied by HT 
series methods. We recall that most existing results on the luiiAtersal combinations of the critical and the correction 
amplitudes are reviewed and thoroughly discussed in Refs.tll€j'c3 

The paper is organized as follows: In section 2 we present our notation and define the quantities we shall study. 
In section 3 we briefly discuss the simplified doubly biased differential approximants which we have used for our 
estimates beside more traditional numerical tools. Our analysis of the series is presented in section 4 along with a 
comparison to earlier series work, to measures perfornaejd in stochastic simulations-ajid-to RG estimates, both by the 
FD perturbative technique and by the Fisher- WilsorcJ e-expansion approacl£i T2j'c3c3. Let ps mention that, very 
recently, the e-expansion of griN) has been extended by Pelissetto and Vicaria from order e^cB to order e^, so that 
we are able to compare our HT results also with their estimates. 

Our conclusions are briefly summarized in section 5. In the appendix we have reported the HT series coefficients of 
X4:{N,P) expressed in closed form as functions of the spixL-dimensionality N. For convenience of the reader, we have 
also reported their evaluation for N = [the SAW modelcH], = 1 [the Ising spin 1/2 model], N = 2 [the XY model] 
andJV = 3 [the classical Heisenberg model]. The present tabulation supersedes and extends the one to order (3^^ in 
Ref.EI which, unfortunately, contains a few misprints. 



II. DEFINITIONS AND NOTATION 



We list here our definitions and notation. As the Hamiltonian H of the A- vector model we take: 

H{v}^~^ E «(^)-^(^')- (1) 

where v{x) is a A-component classical spin of unit length at the lattice site with position vector x, and the sum 
extends to all nearest neighbor pairs of sites. 
The susceptibility is defined by 

x(A,/3)-^M0).i;(f)), (2) 

X 

where (w(0) • v{x))c is the connected correlation function between the spin at the origin and the spin at the site x. 
If we introduce the reduced inverse temperature (A) — 1 — (3 / (if (A) , (here and in what follows 4h stands for 
either sc or bcc, as appropriate, then x{^^ P) is expected to behave like 

X*{N,(3) ~ C#(A)(r#(A))-^W (l + a#(A)(r#(A))^(^) + ... + e#(A)r#(A) + ...) (3) 

when T'^(A) J, 0. C^{N) is the critical amplitude of the susceptibility, a^(A) is the amplitude of the leading 
singular correction to scaling, 0{N) is the exponent of this correction (also called confluent singularity exponent) 
and e^(A) is the amplitude of the leading regular correction. The dots represent higher order singular or analytic 
correction terms. The confluent terms result from the irrelevant variablesE^. Let us recall that not only the critical 
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exponent ^{N), but also the leading confluent correction exponent 0{N) is universal (for each N). On the other hand, 
the critical amplitudes C^{N), a^(A^), e^{N), etc. are expected to depend on the parameters of the Hamiltonian and 
on the lattice structure, i.e. they are non-universal. Similar considerations also apply to the other thermodynamic 
quantities listed below, which have different critical exponents and different critical amplitudes, but the same leading 
confluent exponent d{N). It is known that 0{N) ~ 0.5 for small values of Nl3 and, in the context of the large 
7V-expansiorO, one can also infer that 0{N) = 1 + 0(1 /N). 

Since we have clearly stated which quantities are universal, from now on we shall generally omit the superscript # in 
order to keep the formulas more legible. Notice also that, since there is no chance of confusion, we have systematically 
omitted the superscript + usually adopted for the amplitudes which characterize the high temperature side of the 
critical point. 

The second moment of the correlation function is defined by 

/i2(^,/3) = ^x2(«(0)-^;(f))e (4) 

In the vicinity of the critical point ^2 is expected to behave as 

^i2{N, P) ~ C^(7V)r-^W-2KJV) (^i + a4iV)r''W + ... + e^{N)T + ...) (5) 

as T I 0. |-^ 

In terms of x ^^id /12, the second moment correlation length ^ is definedca by 

In the vicinity of the critical point ^ is expected to behave as 

e(iV, P) ~ Q(iV)T-'^(^) (1 + a4(Ar)r^W + ... + e^{N)T + (7) 

as T I 0. 

The second field derivative of the susceptibility is defined by 

x,y.z r—1 

Notice that this definition differs by a factor from that used in RefB. 

It is well knownE£lEj that, for ^ cx3 at fixed f3 = /3/N, x{N,(3) has a finite non trivial limit x(/3). On the other 
hand, as expected, in the same limit we have xa{N,(3) = 0{1/N). It is the quantity Nx4:{N,f3) that has a finite limit 
Xi{P) simply expressed as 

X4(/3) = -6x^(/3)(x(^)+/3^). (9) 

dp 

Also the — > limit, at fixed /3, existsil and the quantity 

UP) = hm X4(^,/3) - -3 5] CAr,ArJ^i+^^ (10) 

^ N1.N2 

has the following interpretatioiHi^: CN1N2 is the number of pairs (w'-^-', w^^^) of self avoiding walks such that uj^^'^ is 
a A'^i-step walk starting at the origin and w^^-* is a 7V2-step walk starting anywhere and crossing w'^'. 
In the vicinity of the critical point X4(iV, /?) is expected to behave as 

X4{N, fi) ~ C'4(^)r-^W-2^-(^) (1 + a4(^)r^^^) + ... + ei{N)T + ...) (11) 

as T I 0. 

In terms of x, C ^-i^d X4 the " dimensionless renormalizcd four point coupling constant" gr{N) is defined as the value 

of 
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at the critical point (3ciN). Here f{N) = is a normalization factor chosen in order to match the usual field 
theoretic definition of gr{NjEi and v denotes the volume per lattice site {v = 1 for the sc lattice and v — 4/3a/3 for 
the bcc lattice). 

In the vicinity of the critical point g{N,P) is expected to behave as 

g{N, (3) ~ 5^(7V)r^W+3''W-2'^^(^) (l + a<,(iV)T'^W + ... + eg{N)T + (13) 

as T I 0, with 

C|(7V)C2(iV) 

The Gunton-Buckinea mM inequality 

3i'(iV) +7(A^) - 2A4(iV) > (15) 

together with the LebowitB inequality X4(iV, /3) < 0, implies that g{N,f3) is a bounded non-negative quantity as 
r I 0. The vanishing of a(^, Pc) is a sufficient condition for Gaussian behavior at criticality, or, in lattice field theory 
language, for " triviality" C3 of the continuum field theory defined by the A^— vector lattice model in the critical limit. 
If Xii^jf^) is nonvanishing and the above inequality holds as an equality (the hyperscaling relation) 

3v{N) +j{N) -2A4{N) ^0 (16) 

then 

g{N, (3) ~ gr{N) (l + a,(iV)r''W + ... + eg{N)r + ...) (17) 
namely g(N,(3) tends to the nonzero limiting value griN) as r J, 0. 

For cheeking purposes it is useful to recall here the large A'^ limits of the critical amplitudes. They have been 
computedE^I long ago 
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C^%oo) = ^^.^^^ = 0.39228768.. (18) 



with /3f (oo) = 0.2527310098.. and 

Ct^'=(oo) = J- = 0.29974101.. (19) 

647r2(/36-((x)))3 

with Pl'^'ioo) = 0.1741504912.. 

Moreover, we recall that, since in the large N limit 

M2 = g/3x' (20) 

where q is the lattice coordination number, we have C^{oo) — {q(3f {oo)C^ (oo) /6)^/'^ . On the other hand, if we 
denote by Cf{oo) the large TV hmit of NCf{N), by (§) we have Cf (oo) = -12{C*{oo))^ and therefore it follows 
that gf{oo) = 1. 



III. ANALYSIS OF THE SERIES 

As mentioned in the introduction, a variety of careful ajaalyseJ2i~ii of the Ising model HT expansions as well as 
our study of the recently extended iV-vector model serieg3, suggest that the non-analytic confluent corrections to 
the leading critical behavior of the thermodynamic quantities, indicated in the asymptotic formulas (||),(|^), (|^), etc. 
exist and should not in general be neglected in computing numerical estimates of the critical parameters. It has 
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long been observedS~0'0 that these corrections reveal themselves as small apparent violations of both univerSjaJittt 
and hyperscaling in a naive pure power law analysis of the critical behavior. However it is also well knowiH'EHI 
that, unless very long HT series are available, extracting simultaneously estimates for Pc, the exponents and the 
amplitudes of the critical and of the-subleading singularity is a difficult and unstable numerical problem. For this 
task the inhomogcncous DA mcthodL3 of series analysis is generally believed to be more effective than the traditional 
and simpler Fade approximant (PA) method, because, at least in principle, it might be flexible enough to represent 
functions behaving like (j)i{x){x — xo)~'^ + 4'2ix) near a singular point xo, where (j)i{x) is a regular function of x 
and 02 (a^) may contain a (confluent) singularity of strength smaller than uj. Unfortunately, in practice, this is not 
completely true: very long series are needed anyway and/or the procedure should be biased by choosing very carefully 
the structure of the approximaats and by giving proper inputs. We have followed here the latter approach. As in 
some of our previous studieaj'Ell, beside more standard procedures of analysis, we have employed a doubly biased 
prescription which assumes that the confluent exponent 9 and the inverse critical temperature (3c are accurately 
known. This procedure seems to perform reasonably well, even with not very long series. We have taken the values 
of 9{N) as estimated by the FD renprpialization group method. More precisely, for < 4, we have used the values 
suggestedr-bjy Guida and Zinn-JustinO, and for A^ > 4, we have used the six loop estimates recently obtained by A. 
I. SokolovE3 and kindly communicated to us before publication. These values are reported in table 1. We also have 
assumed thatnthe critical temperatures (3f{N) have been determined accurately enough in our previous study of the 
susceptibility^. 

Let us now recall in some detail the features of the simplified DA method. 

We wish to approximate some function, given as a series expansion around f3 = and expected, when /3 "f /3c, to 
have the form 

n=0 

We assume that (3c and the real positive exponent 9 are accurately known, and that 6(/3) and c(/3) are analytic at 
(3 — (3c- We set h{(3c) — bo and c{(3c) = cq. 

We shall estimate the function f{(3) and therefore the constants bo and cq by the particular class of first order 
inhomogeneous differential approximants F{P) defined as the solutions of the equations 

Q™(/3) [(1 - (3/(3c)^^ + ^F{(3)\ + R^{(3) = (22) 
L d(3 (jc i 

with the initial condition F(0) — /q. 

Qm{(3) and Rn{(3) are polynomials of degrees m and n respectively, whose coefficients are calculated, as usual, by 
imposing that the series expansion of F{(3) agrees with that of f{(3) at least through the order in addition 

the normalization condition (5m(0) = 1 is imposed. Assuming for simplicity < < 1, /(/3c) = &o is estimated as 



An.m) _ (3cRn{(3c) 

° ^ --.(/3c) ^"^^ 



and the amplitude co of the subleading term in Eq. (B]]) is estimated by the formula 



where 



^'~Qr,^{t) Q™(/3c) 

We shall consider only the "almost diagonal" approximants with |to — n| < 4. 

The approximants defined by (p2|) are just a small subclass of the general first order inhomogeneous DA's 



(1 - (3/(3c)Qm(fi)^^ + P/(/3)F + i?„(/3) = (26) 

biased with (3c and with 9 by imposing Pi(Pc)/Qmi(3c) — w- Still assuming < < 1, we can estimate &o and cq 
from (PQ) as follows 
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I) _ Rn{f3c) _ PcRniPc) 



(27) 



Jm;n,l) 
-0 



fo 



where 



9 



0c 



g(^-'-''Ht){l-t/(3,y+<^ 
1 , dt 



v^(™;„,0(i) 5(™;".0(/3,)^(l-t//3,) 



i+e 



{(3) — exp 



dt 



(t) QmiPc)' il-t/(3c) 



(28) 



(29) 



are recovered from the general 
= 4-QmiP)- This prescription, 



and D('";"'0(i) has the same form as (p5|). The simple formulas ( ^ ) and (^ 
formulas ( |27| ) and ( |2^ ) by subjecting P;(/9) to the further strong constraint Pi{(3) 
which, for short, we will refer to as simplified differential approxiipaiita (SDA's) might also be viewed as a simple 
DA-like generalization of the biased PA method introduced in Refs.E3~Eil. 

We have carried out many numerical experiments on simple model series having the analytic structure (pT]). They 
show that the SDA's, when biased with the exact values of /3c and 0, are able to produce very accurate estimates of 6o 
and fairly accurate estimates of the confluent amplitude cq. In practice however, we do not have strict control on the 
series: only approximate values of /3c and 9 are available for biasing the SDA's and we do not know the strength of the 
subleading correction terms and of the smooth background. Therefore it is important to understand how sensitive are 
the estimates of &o and cq to the errors in the biased inputs and how they depend on the structure of the singularity. 
It turns out that the estimates of bo are rather stable when the biased value for j3c and for 9 are varied away from 
their true values in a range comparable to the typical estimated uncertainties in the realistic cases. On the other 
hand, cq appears to be much more sensitive to errors in the biased values. Let us consider, to be definite, the case of 
the very simple test series 



/(/?) 



Pc 



boexp{l - ^) 

Pc 



+om-(3f) 



(30) 



which we have examined for various values of 9. If the size of the subleading correction to scaling is much smaller than 
the size of the leading one, namely if | ci |<c| cq | and we bias the calculation with the exact values of the parameters 9 
and /3c, we are able to determine bo by (p3| ) to within less than 10~^% and cq by (|2^ ) to within less than 1%. However, 
if the SDA's are biased with a value of 6 which is off the right value by 5%, then the relative error of cq can become 
as large as 15%, while the error of bo increases to some 0.1%. The precision of bo remains essentially unchanged, but 
the sensitivity of cq to variations in the biased values and, as a consequence, the accuracy of its estimate is somewhat 
worsened in the slightly more complicated, but sometimes realistic case in which | ci |«| co |. Unsurprisingly, the worst 
situation occurs when the leading confluent amplitude is much smaller than the subleading one, since the uncertainty 
in the numerical estimate of cq may then become very large. In conclusion, taking a conservative attitude, we can 
safely expect that, for the HT series we are going to study, the relative error on the value of /(/3) at /3c can be much 
smaller than 1%, while the uncertainty of the correction amplitude can be as large as 20%, unless the amplitude is 
very small: in this case, due to a higher sensitivity to the biased values and/or to the neglect of possibly important 
subleading corrections, our estimates are likely to be much more inaccurate. In order to better understand these 
results let us also observe that, if we tried to estimate bo in (^) by simple PA's biased with /3c, the relative error 
would be substantially larger and increasing with the size of the correction amplitude. Finally, we remark that in all 
computations presented below, the error estimates are always somewhat subjective. They include effects both from 
the scatter of the approximant values, possible residual trends in sequences of estimates, as well uncertainties of the 
bias inputs. 

We have applied the SDA approximation procedure not only to the quantity g{N,j3) in order to compute the 
confluent amplitude ag{N), but also to the "effective exponent" of Xi 



74(iV, /3) ^ - = ^IN) + 2A4(iV) - a4N)9iNy^^^ + o(/W) 



(31) 



in order to compute the critical exponent and the confluent amplitude ai{N). 
Moreover we have examined the analogous quantities 
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7(iV, P) ^ (/3,(iV) - /3)^^^^MpZZ ^ ^(iv) - a,(iV)^?(7V)r^W + o(r«W) (32) 
in order to compute the confluent amplitude a^{N), and 

.(/3, A^) ^ - ^ ,(^,) _ adN)e{Ny(^) + o{rO(^^) (33) 

in order to compute a^{N). Notice that the estimates thus obtained for the confluent amplitudes a^, a^, and 04 
are biased solely with f3c and 6. However, due to their definition as residua, the sensitivity of the results to the biased 
value for /3c is higher than in the case of gr- 

The estimates of the critical amplitudes have been obtained by examining quantities like 

T-^Wx(^, P) ^ C^{N) (1 + a^(7V)/(^) + ... + e^{N)T + ...) (34) 

or the analogous expressions for X4 and This procedure also yields the correction amplitudes, but since it requires 
biasing also with the critical exponent ^{N) (or v{N) etc.), we expect that the corresponding results will be subject 
to a larger uncertainty. 

In conclusion, whenever sizable confluent corrections are present, the doubly biased SDA procedure will produce 
values of gr{N) which are slightly, but definitely different from estimates by generic DA's not directly constrained to 
reproduce the confluent singularity and, a fortiori, from the simple PA estimates. Indeed, since 6 < 1, the function 
g{N,(3) will approach with a divergent slope its value at j3c{N), from above if the correction amplitude is positive 
or, otherwise, from below. As a consequence, too smooth extrapolations of g{N,j3) to the critical point Pc would 
overestimate the correct result in the former case and underestimate it in the latter. Analogous problems will occur 
in the study of the exponents and of the correction amplitudes for X; M2j X4) the only difference being that, since in 
the formulas for the effective exponents (^l|),(|3^) and (^) the correction amplitudes appear with a negative sign, the 
critical exponents will be overestimated if the amplitudes are negative and they will be underestimated otherwise. 

Let us add finally that throughout our work we have not relied solely on the above numerical technique, but we also 
have always considered various other approximations obtained by more conventional methods in order to understand, 
or at least to be aware of any differences in the estimates. 



IV. RESULTS AND COMMENTS 



Since our analysis is aimed at exposing the role of the non analytic corrections to scaling, it is desirable firstly to 
test whether the values of the confluent exponents taken from the FD perturbative computations are also generally 
consistent with the estimates, unfortunately not yet as precise, which can be extracted directly from the HT series. 
Indeed, as we have mentioned above, the amplitudes of these corrections are not universal and therefore they might 
be negligibly small. One might even suspect that our analysis is somehow artificially forcing on the series a behavior, 
which, due to their insufhcient length, they are not yet able to exhibit. On the other hand, it has been argued that 
the uncertainties usual ly mioted for the FD values of the renormalized couplings and of the confluent exponents might 
be unrealistically smallEilc3'EJ'E3. In fact, one should recall that in the context of the three-dimensional \{<PY ^^1^ 
theory, the confl ueut ftaponent is computed in terms of the slope of the beta- function at the fixed point gr{N). As 
indicated in Refs.Elliy'C3, the presence of non-analytic terms, with sufficiently large amplitudes, in the expansion of 
the beta- function at gr{N), might spoil the convergence of the estimates both of the renormalized couplings and 
of the confluent exponents. The ensuing uncertainties would reflect on the accuracy of the estimates of the critical 
exponents. Moreover the g— expansion of the critical exponents would itself be directly affected by similar non-analytic 
contributions. The pragmatic point of view adopted in RefEj is that if these singular terms exist, they do not seem 
to have visible effects. 

Let us then show that the values of 0{N) reported in Table 1 are approximately consistent with thcractual behavior 
of the series. Assuming knowledge only of (3f{N), we have computed the Baker-Hunter transforma^ of the x and 
/i2 series and, by reconstructing the locations and the residua of their singularities, we have estimated exponents and 
amplitudes of the critical singularity and of the leading correction to scaling. Unfortunately this procedure fails to 
detect narrow and clear signals of the scaling corrections for < 4, probably due to the small size of their amplitudes. 
But the situation is completely different for > 4. In this range of values of N, the Baker-Hunter method leads 
to values of d{N) fairly consistent with those reported in Table 1. Also the values of the correction amplitudes, are 
compatible with those emerging from the SDA analysis to be discussed below. Moreover, the results are rather stable 
in a relatively wide range of biased values for Pc- We regard this as convincing evidence that the confluent corrections 
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cannot be a by-product of our double biased analysis and as a confirmation that their amplitudes are not small for 
> 4. Unfortunately, the uncertainties which affect this method for estimating the confluent exponents and the 
correction amplitudes are still rather large. For instance, using the bcc lattice series for x, the Hunter-Baker procedure 
suggests 

61(4) = .64(4); 6(6) = .63(4); 61(8) = .66(4); 61(10) = .69(4) (35) 

A second consistency test can also be performed. On both lattices and for each value of N, we have studied how 
our SDA estimates of 5,. (N) depend on the biased value used for the confluent exponent by varying it in a 20 — 30% 
range around the central value 9{N) indicated in Table 1. For all values of N such that the confluent amplitudes are 
not too small, it has been quite interesting to observe that, although the estimates of griN) obtained from the sc and 
the bcc series are in general somewhat different for a generic value of 9, they tend to become equal, or at least very 
close, when 9 ~e{N). 

These two tests give us further confidence that the main lines of this analysis and the specific biased values of 9 
used as inputs are reasonable. 



A. Hyperscaling tests 

We shall now proceed to examine directly X4(iV, P) in order to estimate its critical exponent 7(A'^) -I- 2A4(iV) and to 
compare it with the value 2^{N) + 3z/(7V) it should take if the hyperscaling relation ( |lq ) holds true. On each lattice, 
the analysis has been performed by first-order SDA's of the effective exponent (|3^) doubly biased with 9{N) and with 
the value of f]f (N) obtained in our previous (biased) analysis^ of the susceptibility. 

We have reported in table 2 our estimates for the critical exponents of x^i^j P) obtained by this procedure together 
with the biased values of Pf{N) and the values of 2j{N) + 3h'{N) obtained by the analogous biasing procedure in 
our previous analysistl of x and No significant violation of universality and hyperscaling is observed. Notice that 
no such extensive test of hyperscaling exists so far in the literature. 

Let us quote a few earlier studies of this issue for particular values, of N. In the = case, a study of X4 
based on (|l^) has been performed by a Monte Carlo simulation in RefO. The authors have measured the exponents 
2A4-7 1.7317 ±0.0074 ±0.0074 and 1/ = 0.5745 ± 0.0087 ± 0.0056. The final result is expressed as 3v + -/-2Ai ^ 
—0.0082 ± 0.0027 ± 0.018, the first error being systematic and the second statistical. |— in 

In the iV = 1 case, the tests of the hyperscaling relation (|l6|) are Bipaerous and have a long historyEj'E^I. The validity 
of ( p^ ) for the 3d Ising model had been questioned by G. A. BakeioE3 on the basis of an analysis of 10-12 term series 
for the sc, bcc and fee lattices. A few years later, when B. G. Nickel computed 0{(3^^) series on the bcc lattice for 
X and /i2 in the spin S Ising model, it became clear that rather long-sepfts-|Were necessary to allow for the scaling 
corrections and thus to obtain more satisfactory estimates of 7 and t^lhE^Pn. .pji the other hand accurate analyses 
of the critical behavior of the X4(l? /?) series to order on the sc latticeEJIljOEZl had^ielded reliable values also for 
A4(l). On the basis of these results, as well as of various recent Monte Carlo resultsEa"L2l a common consensus was 
reached that, for = 1, if any violation of (|6|) occurs, it should be much smaller than was initially suspected. Our 
contribution to this issue also consists in providing an extension from order f3^^ to order /3^^ of the Ising bcc series 
for X4, and therefore in further improving the accuracy of the HT test of hyperscaling and universality even for the 
widely studied = 1 case. 



B. Renormalized couplings 

Let us first mention that, since — 0(/3) in the vicinity of P — 0, from the series for x, and X4 we can form 
two distinct auxiliary functions w{N,P) and u{N,P), analytic at /3 = 0, both of which, when extrapolated at /3c yield 
g{N,Pc) and therefore gr{N), if we assume the validity of the hyperscaling relation. More precisely we shall consider: 

(./(A^)x4(iV,/3))2/3 ^^^^ 

whose value at Pc{N) is ^^(A^)^^/'^, and 

W(N,P) = :r-p, ■ (37) 
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whose value at (3c{N) is gr{N). 

It is interesting to form approximants both of u(/3) and of u'(/3) because for various values of N, at the presently 
available order of expansion, they still show slightly different convergence prpperties. This may be seen as an indication 
that the Xi series are still not very long. Indeed, as we have argued in Ref.Ell, at order the dominant contributions 
to the HT expansion of come from correlation functions of spins whose average distance is ~ s/4. Therefore the 
presently available expansions with Smax = 17 still describe only a rather small system. 

Table 3 contains our estimates of the universal renormalized coupling gr{N). 

For TV < 4 we have evaluated gr{N) by forming SDA's of the auxiliary function w{N,P), which has been chosen 
because it yields sequences of estimates showing little or no residual trends when an increasing number of series 
coefficients is used. On the other hand, for > 4, we have used u{N, f}) because the estimates obtained from it 
show the slowest (generally decreasing) residual trends. Whenever relevant, we have indicated this fact by reporting 
asymmetric error bars. 

In the = case, allowance for the correctiopr,to scaling yields a value of 5r(0) approximately 2% smaller than the 
one recently obtained within the FD expaaaioriEl, but very close to the value suggestedca by the e— expansion. Our 
value is also close to that indicated in Ref.c2l and produces, via the seven loop FD perturbation series, central values 
of 7(0) and viQ) « 0.2% lower than those quoted in Ref.Ell, but within their error bars. It is also worth recalling that 
also our earlier HT analysis^ of x(0,/3) and ^^(0,/3) had supported those low exponent estimates in good agreement 
with very recent high precision measures by stochastic methods on the sc latticeO. 

For iV = 1, on the sc lattice, we have reported here a central estimatepof 5, (1) slightly lower than, though consistent 
with the estimate (?r(l) = 1.411 obtained from our previous analysisEll based on SDA's of u(l,/3), rather than of 
«;(!,/?). 

A small sample of the most recent estimates of 5r(l) by various methods has also been included in the table. All 
of them appear to be mutJialily consistent, if we consider how difficult it has been to achieve very accurate, Monte 
Carlo measures of (7^'^(1)E3"l2I and we recall that, even in the Ising case, the previous HT series estimates^ of the 
renormalized coupfing were based on expansions shorter than those presented here. Indeed, although X4(l,/3) on the 
sc lattice has long been known thmugh order the correspapding expansion for the renormalized coupling was not 
available before our recent workoED, because ^^(1, (3) reachedlla only order (3^^ . On the other hand the X4(1j P) series 
for the bcc lattice was known to order (3^^ only. 

To our knowledge, no Monte Carlo results are yet available for A^ > 1. [—in 

For A^ > 3 our estimates are systematically slightly higher than the FD values of RefsHH, and perhaps the 
residual decreasing trend in our estimates might not be sufficient to reconcile them. This difference is related to our 
allowance of the scaling corrections byi-doubly biased SDA's and is consistent with the higher values of 7 and v that 
we had obtained in our biased analysisB of x and f . As we have stated above in discussing the general features of the 
SDA's, significantly larger estimates for A^ < 4 and somewhat lower estimates for A^ > 4 would be obtained, if the 
renormalized couplings were evaluated by simple PA's. This fact is completely consistent with the observed behavior 
of the correction amplitudes as functions of A^ to be discussed in next s.ubsection. A similar observation has been made 
also in Ref.o where, on the basis of the old sc lattice 0(/3^'*)|-seEiea3^Q, the gr{N) have been evaluated by ordinary 
DA's, either directly or after performing a change of variableO^Ell designed to regularise the leading correction to 
scaling and numerically similar to our SDA's. Therefore the final HT estimates of Ref.O essentially agree with ours. 

We. have included in table 3 some estimates of gr{N) based on the e-expansion to order recently presented in 
Ref.O. They are compatible with ours for A^ < 3, while, for A^ > 3, the central values are ~ 2% lower. 

C. Critical and correction amplitudes 

In tables 4 and 5 we have reported our estimates of the (non-universal) critical amplitudes C'^{N), C'^'^{N) and 
C^^iN), Cf^iN), based on the values of I3*{N), -f*{N) and v*{N) obtained in the biased analysis of Ref.i. Earlier 
determinations of the critical amplitudes either from the extcapolation of (gcnerallj!-, shorter) HT series or from 
stochastic simulations are also available for A' = 1 in Refs.OO, for A^ = 2 in Ref.LJ, and for A^ = 3 in Ref.O. 
However, comparisons with the results of tables 4 and 5, which in general are close to the earlier ones, are not very 
illuminating, since the estimates depend sensitively on the numerical procedures, on the biased values used for (3f{N) 
and on the relevant critical expouejats, which are slightly different in the various studies. 

For example in the A^ = 1 caseLZrE3 (on the basis of shorter sc lattice series, but the same bcc series), the following 
estimates are proposed: C^'=(l) = 1.0928(10) and Cf (1) = 0.4984[5°|; C^^^il) = 1.0216(8) and ^=(1) = 0.4608(2) 

assuming /3f (1) = 0.221630(12), /3^'='=(1) = 0.157368(7), 7(1) = 1.2395 and = 0.632. 

For Af = 2, in Ref.Lfl the estimates C^^{2) = 1.058(7) and C|'=(2) = 0.498(2) have been obtained from a fit to 
Monte Carlo data, assuming 7(2) = 1.3160(25), v{2) — 0.669(2), allowing for confluent corrections with exponent 
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e{2) = 0.522. This fit also yields /3f (2) = 0.454162(9) from the analysis of x, and /3f (2) = 0.454167(10), from the 
analysis of ^. |— ■ 

For 3, in Ref.O, the estimates CJ'=(3) = 0.955(6) and C|'=(3) = 0.484(2) have been obtained from a fit to 

Monte Carlo data (with no allowance for confluent corrections) also yielding = 0.69294(3), 7(3) = 1.391(3) 

from the analysis of x and /3f (3) = 0.6928|li4), t/(3) = 0.698(2), from the analysis of C 

As has been stressed in general in Ref.t3 and as we have anticipated in our considerations of Section 3 on the 
numerical properties of SDA's, the discussion of the estimates of the scaling correction amplitudes is much more 
delicate. Let us first comment on some qualitative features of the estimates of these amplitudes for the sc and the 
bcc lattices which are denoted as a'^%N), a^^^^N), af{N), af''{N), a''g'={N), a^g'"'{N), and reported in table 6. 

Both correction amplitudes a^{N) and d^{N) are negative for N ^ 2, whereas they are positive and increasing 

for N > 2. (Actually we have reported a positive value for a'^'^{2), but with a large uncertainty.) Therefore the ratio 

a^/ux is very likely to be positive for all values of N. The correction amplitudes a'^{N) and a^{N) turn out to be 
small, but not negligible for < 1 and rather large for iV > 4, both in the sc and in the bcc lattice case. On the 
contrary, for N = 2 and = 3 they are very small. Thus the overall behavior of the correction amplitudes for 
and ^ as functions of N appears to be smooth and completely consistent with the size and the sign of the differences^ 
between our unbiased estimates of the critical exponents 7 and v and the corresponding estimates biased with both 
Pc and 9. More precisely, we recall that the non-analytic corrections to scaling lead to (slightly) higher effective 
exponents for < 2 and to (significantly) lower effective exponents for A^ > 4. On the other hand, the a^(A^)'s are 
positive for N ^ 2, while they are negative and decreasing for A^ > 2, so that the ratio ag/a^ is negative for any A'^. 
(Actually we have reported a positive value for ag^{'i), but with a large uncertainty.) The af{Nys are generally not 
small, except for A^ = 3 in the sc case and for A^ = 2 in the bcc case. 

It is appropriate now to quote some eajlie.r evaluations of and by HT series or Monte Carlo simulations. In 
the TV = 1 case, it had beenjestablishedEj^L^HO long ago that the sign of ax{l) and of a^{l) is negative on the sc, 
bcc and fee lattices. In Ref.LJ, for the spin 1/2 Ising model on the bcc lattice, the estimates a'^'^{l) w —0.13 and 

a|'^'^(l) ~ —0.11 have been indicated, together with the central values 7(1) = 1.237, = 0.630 and ^1) = 0.52, on 
the basis of a second order DA analysis. For N — 2, the above cited Monte Carlo simulation of Ref.O, yielded the 
estimates a^i2) = —0.15(6) and a|^(2) — —0.20(4). Clearly, in both cases the critical parameters are slightly different 
from ours and this is sufficient to explain the somewhat different estimates for the correction amj^iitudes. 

In the spin 1/2 Ising case, it has been argued long agotil that ag^^(l) should be large. Recently^, it also has been 

observed that, if the sc lattice Monte Carlo data of Ref0 are simply fitted by the function griP) = 3*(1 + o-gi'^^'^), 
the value ag'^(l) ~ 1.13 is obtained in fair agreement with our own estimate. 

In table 7 and 8 we have listed some earlier estimatfts-|pf.|thje-. universal ratios a^{N)/ax{N) and ag{N)/ax{N) 



of correction amplitudes obtained by various methodscJacaSLZTLa. We believe that, for A^ < 4, it is not very 
meaningful to quote the ratios of our central estimates of a^, and ag. Indeed, as we have already pointed out, for 
these values of A^, the amplitudes a^, are small and very sensitive to the biased inputs. As a consequence, these 
parameters must be finely tuned, which cannot be justified until longer series will be computed. We shall indicate 
below a possible alternative way out o£,this difficulty. However, the case A^ = 1 deserves further commeaL In this 
case, on the sc lattice, a very accurateE3 determination of /3c^(l) is available, and also the value of /3J^^(1)l3 appears 
to be sufficiently safe, so that the ratios of our central estimates of the amplitudes are more trustworthy and we have 
reported them in parentheses. 

It is interesting to recall also that, for A^ = 1, suggestions that ag/a^ should be large came both from earlier HT 
estimatesLj on the fee lattice (a^'^'^(l)/a^'^'^(l) « 3.9) and from the RG estimates reported in Table 8. This is a further 
hint that the corrections to scaling should not be neglected in computing 5r-(l)- r-in 

Unfortunately, the e-expansions of these universal ratios presently only extend to second orderQ'ES, so that again 
we have to point out that the uncertainty of the corresponding estimates might be larger than, indicated. As we have 
mentioned above, even the estimates of these ratios from the much longer FD expansionsEEEj might have problems. 
For A^ < 4, as already observed, all series including ours are too short to accurately extract the correction amplitudes. 
This is particularly the case for Xi- Moreover, when longer series become available, our approximation procedures 
might need some improvement. Nevertheless these first results from HT series on an extended range of values of N 
seem to be qualitatively very reasonable. 

As an indication of work in progress, we wish to add that, even within the present order of expansion, somewhat 
more accurate estimates of the critical parameters, are likely to be obtained by proceeding systematically in the spirit 
of the Chen, Fisher, Nickel and Rehr approachE3E3. In the A^ = 1 case on the bcc lattice, these authors have examined 
HT series for families of models specified by an appropriate continuous auxiliary parameter. The members of these 
families interpolate between the spin 1/2 Ising and the Gaussian model and all of them are good candidates for 
belonging to the same universality class. (This approach easily generalizes, in various ways, to other values of A^ and 
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it is a virtue of the LCE method that the corresponding series can be derived essentially with no further computational 
effort.) By varying the auxiliary parameter, these authors have selected representative models such that the leading 
correction amplitudes a^, (and ag) vanish. It is clear that, under these conditions, even by employing ordinary 
unbiased DA's, the accuracy of the estimates of the universal quantities can be improved dramatically. On the other 
hand, within the same approach, it is probable that also the correction amplitudes will be more accurately measured 
by focusing on representative models in which they are sufficiently large, provided, of course, that the subleading 
terms are not even larger. Thus more reliable estimates might be achieved for their universal ratios, in particular in 
the range TV < 4. 

V. CONCLUSIONS 

The main result of this paper is the extension through 0{(3'^'^) of the series for X4(iV, for arbitrary N, on the 
sc and on the bcc lattices. Both sets of expansion coefficients have been tabulated in the appendix in order to make 
independent checks of their correctness and alternative analyses conveniently feasible. 

A second interesting result is the numerical analysis of the critical behavior of X4(A^)/3) which confirms fairly well 
the validity of universality and hyperscaling over a wide range of values of N. We have also presented a first estimate 
of the size of the scaling corrections for x, and X4 and, allowing for them, we have improved the accuracy in the 
determination of the critical amplitudes and of the renormalized couplings. 

The agreement between our estimates of gr{N) and those from the RG approaches is generally fair, but not always 
perfect. At this level of approximation, it is premature to emphasize such minor discrepancies. We believe, however, 
that longer HT series for all quantities studied here and perhaps improved analyses are still of some interest to achieve 
more reliable estimates and to reduce the error bars substantially. Considering the performance of our codes, these 
are presently quite realistic objectives and, therefore, work is presently in progress to compute further expansion 
coefficients. 
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TABLE I. Values for 8 used in our biased evaluations and determined by FD perturbative expansion 



N 





1 


2 




3 


4 


6 


8 


10 




9 


.478(10) 


.504(8) 


.529(8) 




.553(12) 


.573(20) 


.626(10) 


.670(10) 


.707(10) 






TABLE IL Verification of hyperscaling 


for various 


values of N. 








N 




Lattice 










7 + 2A 




27 + 3;^ 


5 







sc 
bcc 




.213493(3) 
.153128(3) r-. 




4.10(2) 
4.081(8) 




4.0822(34) 
4.0801(34) 


1 




sc 
bcc 




.2216544(3)y 
.157373(2)E 




4.361(8) 
4.366(6) 




4.3721(44) 
4.3692(27) 


2 




sc 
bcc 




.45419(3) 
.320427(3) 




4.665(20) 
4.663(15) 




4.675(12) 
4.666(12) 


3 




sc 
bcc 




.69305(4) 
.486820(4) 




4.953(20) 
4.948(15) 




4.960(12) 
4.946(12) 


4 




sc 
bcc 




.93600(4) 
.65542(3) 




5.24(2) 
5.22(2) 




5.259(17) 
5.236(17) 


6 




sc 
bcc 




1.42895(6) 
.99644(4) 




5.67(2) 
5.65(2) 




5.691(19) 
5.673(17) 
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TABLE III. The renormalized coupling constant gr{N) for a range of values of A*" on the sc and the bcc lattice as obtained 
by various methods. 



iV 



HT sc 



HT bcc 



e-exp. 



1.390(17)&" 
1.397(8)0 
1.413(13)0 
1.387(7)0 
1.366(15)0 



FD ex 




Monte Carlo 





1 

1 
2 
2 
3 
3 
4 
4 
6 
8 
10 



1.388(5) 

1.408(7)p 
1.459(9)0 
1.411(8) 

1.409(10) 

1.392(10) 



J--300/ 10) 

1.320 +^>, 



1.387(5) 
1.407(6) 
1.411(6) 
1.406(8) 
1.394(10) 

1.326 +f,) 
1.294|+fg) 



1.295(7) 







1.391(30) i 
1.462(12)Ea 



TABLE IV. Critical amplitudes on the sc lattice for various values of N. 



N /3c (A) j{N) iy{N) C°"(A) Cf (A) 

0.213493(3) 1.1594(8) 0.5878(6) 1.115(1) 0.5101(3) 

1 0.2216544(3) 1.2388(10) 0.6315(8) 1.111(1) 0.5027(3) 

2 0.45419(3) 1.325(3) 0.675(2) 1.014(1) 0.4814(3) 

3 0.69305(4) 1.406(3) 0.716(2) 0.9030(8) 0.4541(2) 

4 0.93600(4) 1.491(4) 0.759(3) 0.7571(8) 0.4155(2) 
6 1.42895(6) 1.614(5) 0.821(3) 0.6054(8) p 0.3708(2), 
oo 2. 1. 0.392287..N 0.314870.. I 



TABLE V. Critical amplitudes on the bcc lattice for various values of A. 



N PcjN) 7(A) liN) C^'^jN) Cf"(A) 

0.153128(3) 1.1582(8) 0.5879(6) 1.087(1) 0.4846(2) 

1 0.157373(2) 1.2384(6) 0.6308(5) 1.034(1) 0.4659(2) 

2 0.320427(3) 1.322(3) 0.674(2) 0.918(1) 0.4371(2) 

3 0.486820(4) 1.402(3) 0.714(2) 0.794(1) 0.4072(2) 

4 0.65542(3) 1.484(4) 0.756(3) 0.6580(8) 0.3691(2) 
6 0.99644(4) 1.608(4) 0.819(3) 0.5020(6) p 0.3231(2), 
oo 2. 1. 0. 299741. 0.263818.. I 



TABLE VI. Correction amplitudes on the sc and the bcc lattices for various values of A. 



N 


a- (A) 




4^{N) 


af'=(A) 


a- (A) 


a^-(A) 





-0.022(10) 


-0.05(3) 


-0.11(3) 


-0.1 


1.5(3) 


1.9(4) 


1 


-0.10(3) 


-0.08(3) 


-0.12(3) 


-0.08(3) 


1.0(2) 


1.0(2) 


2 


-0.04(2) 


0.01(2) 


-0.07(3) 


-0.005(9) 


0.39(8) 


0.14(3) 


3 


0.06(3) 


0.17(3) 


0.003(6) 


0.09(3) 


0.05(10) 


-0.29(6) 


4 


0.30(6) 


0.5(1) 


0.14(3) 


0.28(6) 


-0.12(3) 


-0.55(10) 


6 


0.73(15) 


1.1(2) 


0.37(8) 


0.60(15) 


-0.35(8) 


-0.77(20) 


8 


1.1(2) 


1.8(3) 


0.53(10) 


0.92(20) 


-0.43(10) 


-1.0(2) 
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TABLE VII. The universal ratios of correction amplitudes a^{N) / a^{N) for various values of A'^. 



iV HT so HT bcc e-exp. FD exp. r-. 

0885(50)0 

1 (l-2(4)) (^-^^^l)] 0.762(3000 

1 O.7O(3)@;O.85(5)0 0-56^5)0 O.65(5)0p 

2 O.eSLlI 0.687(10*3 

2 0.615(5)p 

3 0.63M5pl 

3 o.6olia 

4 0.49(15) 0.55(15) 
6 0.51(15) 0.55(15) 
8 0.49(15) 0.54(15) 



TABLE VIII. The universal ratios of correction amplitudes ag{N) / a-^[N) for various values of N . 



N 



HT sc 



HT bcc 



e-cxp. 



FD exp. 



3.0(5) 



-1.2(4) 
-1.1(4) 



-3 



5(5)) 



-2.2(5)1 







-2.85(6) 

-2.08(5)2 
-1.65(4)r 



-0.7(2) 
0.8(3) 
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APPENDIX A: THE SECOND FIELD DERIVATIVE OF THE SUSCEPTIBILITY ON THE SC LATTICE 



The HT expansion coefficients of the second field derivative of the susceptibility 
on the sc lattice are 



di{N) 



di{N) 



di{N) = -48/Af 

, , , -1248 - 660 AT 
= 7V3 (2 + TV) 

, , , -12480 - 6912 AT 
'^'^''^ = jV4 (2 + iV) 

-851712 - 1128192 N - 474000 - 61236 JV^* 
JV5 (2 + N)^ (4 + N) 

-6573312 - 8786880 JV - 3725856 JV^ - 483840 
iV6 (2 + JV)^ (4 + JV) 



-565908480 - 1137490944 iV - 877991616 iV^ - 321566208 N'' - 55124016 iV* - 3514968 iV'* 
^ ' JV7 (2 + Nf (4 + AT) (6 + N) 

For the coefficients which follow it is typographically more convenient to set dr{N) = Pr{N)/Qr{N) and to tabulate 
separately the numerator polynomial Pr{N) and the denominator polynomial Qr{N), 

Pt{N) = -3849744384 - 7739114496 iV - 5976661248 AT^ - 2190260352 AT^ - 375495552 AT* - 23938560 AT^ 

Qr{N) = AT* (2 + Nf (4 + AT) (6 + AT) 
Pa{N) = -1607361822720 - 4630934495232 AT - 5658731108352 - 3815032300032 - 1545703906176 AT* 



QaiN) = at'' (2 + AT)* (4 + Nf (6 + AT) (8 + AT) 
Pg{N) = -10146634334208 - 29145299066880 Af - 35515117682688 AT^ - 23883563143680 JV^ - 9654774400512 AT" 
-2393329321728 AT^ - 354292592640 AT** - 28580124768 AT^ - 960719616 AT* 
Qg(N) = (2 + AT)" (4 + Nf (6 + AT) (8 + AT) 
Pio(Ar) = -4986778413957120 - 18502905604472832 AT - 30434129019666432 - 29229711376023552 
-18173047179460608 Ar"* - 7663189901412864 AT** - 2231748901824768 AT** - 448039233434880 N'^ 
-60661040264064 AT^ - 5267106682272 - 263609235360 AT ^° - 5754914568 AT", 
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Qio(JV) = (2 + N)" (4 + Ny (6 + AT) (8 + N) (10 + JV) 
Pii{N) = -29957292231229440 - 110697821461807104 AT - 181359443339575296 - 173526582721855488 Af^ 
-107505154356572160 AT* - 45183600633858048 - 13119088167641088 Af" - 2626537796155392 N'' 
-354741723247872 AT^ - 30735239578752 Af^ - 1535369868288 Af^° - 33465664512 Af^"^ 
Qii (N) = N'^^ (2 + Nf (4 + Nf (6 + AT) (8 + AT) (10 + N) 
(N) = -25424963775485706240 - 112568160566969892864 N - 226122593284806672384 - 272739288108751650816 AT^ 
-220335024310830366720 AT* - 125926002861175824384 AT^ - 52429589825842464768 AT** - 16133360608103531520 n'^ 
-3682631735427354624 AT** - 619886922488017920 Af" - 75677317494258048 Af^° - 6492816409650048 Af" 
-369850295426784 Af^^ - 12514445149200 Af^^ - 189708636600 Af" 
Qi2(Af) = Af^^ (2 + Nf (4 + Af)^ (6 + Nf (8 + N) (10 + Af) (12 + Af) 
Pi3{N) = -147436255220890337280 - 649808498848293715968 N - 1299512129528339103744 AT^ 
-1560658843095801004032 AT^ - 1255571803469692796928 A?* - 714760376636548620288 N^ 
-296490805364682670080 Af" - 90920656735938183168 n'' - 20688048472922093568 AT* 
-3472336842523894272 Af" - 422813491673485824 Af^" - 36192084771724800 At" 
-2057419542670080 Af^^ - 69492377025984 Af^^ - 1051829121024 AT^* 
Qi3(Af) = Af^'' (2 + Nf (4 + Af)^ (6 + Nf (8 + N) (10 + Af) (12 + Af) 
Pii{N) = -566463587862368653148160 - 3043875242905724409348096 N - 7576777562322913437155328 
-11598723404301679271608320 - 12226130280460910006894592 Af* - 9415725334597698602139648 Af^ 
-5486072371880331913592832 N" - 2470783134534602565992448 Af'' - 871335330057590064906240 AT® 
-242243678099334349271040 Af" - 53184931154543324258304 Af^" - 9193953185463740998656 Af" 

-1241394935192535991296 Af^^ - 129075746454058556160 AT^^* - 10102244514482628864 Af" 
-573981185124034560 Af'^'^ - 22283021804865984 Af^** - 527221923353952 Af^^ - 5719330613520 Af"^* 
Qii{N) = N^^ (2 + AT)'' (4 + AT)* (6 + Nf (8 + N) (10 + AT) (12 + N) (14 + AT) 
Pi5(Af) = -3199047998804219563868160 - 17113814040741606825394176 Af - 42413540836361700246552576 N^ 
-64650199580585358413266944 AT^ - 67863921600679932441133056 AT* - 52054127611661734870253568 AT^ 
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-30212418901698609794777088 JV° - 13556920252810313722822656 JV'' - 4764324172851484585525248 AT* 
-1320233119260810969808896 - 288978328415607033544704 - 49814490389339231121408 JV" 
-6708703760313406396416 JV^^ - 695898456216849782784 JV^^ - 54348379646212915200 JV" 
-3081948095561442816 JV^^ - 119439454925014272 JV^** - 2821599683684352 JV^'' - 30566943406080 AT^® 
Qi5(JV) = JV^^ (2 + Nf (4 + AT)* (6 + Nf (8 + JV) (10 + N) (12 + JV) (14 + N) 
Pie{N) = -18299002328234221207226941440 - 114585264379209072949397028864 N - 336677952680726802363970486272 
-617145013898356017378109685760 JV^ - 791521185299488842160294330368 AT* 
-755125087023410853999229796352 n'^ - 556163233950005337533492232192 N*^ 
-323972815609871396650371514368 AT'' - 151693834060073149464108859392 
-57712428469432137886012145664 - 17963030647495383412208861184 Ar^° 
-4590796640925368216011948032 AT" - 964233613788142396008136704 AT^^ 
-166150396381094082218360832 Af" - 23381941196613066788198400 
-2666695999548097875723264 AT^^ - 243616432556930571577344 AT^** 
-17525072840959462713984 Af^'' - 968152139098787129664 AT^* - 39538186077707686464 AT^^ 
-1121781497993998176 Ar^° - 19697772998733576 AT^^ - 160868281485204 AT^^ 
Qie(Ar) = AT^^ (2 + Nf (4 + Nf (6 + Nf (8 + Nf (10 + N) (12 + N) (14 + AT) (16 + N) 
Pij{N) = -101235958788462673253105664000 - 631215419419149274713371443200 AT - 1846824461720939657841963171840 AT^ 
-3371230360792179922631399571456 N^ - 4306156546428529311715461955584 AT* 
-4091812122737075444969361113088 AT^ - 3002060652157587317256634761216 AT" 
-1742224189122069022527046287360 AT^ - 812841546356917149508011294720 N^ 
-308190677384300600694770368512 N^ - 95612813294752704503595663360 N^" 
-24360624222950336197451317248 A?" - 5101833629477951341726236672 N^^ 
-876738435402445320757764096 N^^ - 123070791802592404991434752 Af" 
-14003418390237375888211968 A?^*^ - 1276532972320935406141440 A?^'' 
-91648782429013822778880 N" - 5053887369689137387008 N^^ - 206055475899900727296 AT^^ 
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-5837513893110173568 Ar^° - 102365100472039776 iV^^ - 834986039880960 JV^^ 

QieiN) = Af^** (2 + Nf (4 + Nf (6 + Nf (8 + N)^ (10 + N) (12 + N) (14 + N) (16 + N) 

In particular for = we have (in terms of the variable j3 = (3/N) 

X4(/3) = -3 - 72/3 - 936;3^ - 9360^9^ - 79848/3* - 616248/3^ - 4421160/3** - 30076128/3'' - 196211160;9* 
-1238602824^3^ - 7609219992/3^° - 45711200304/3" - 269412610536/3^^ - 1562290776792/3^^ - 8932238341992/3^* 
-50443946980992/3^^ - 281783630311272/3^" - 1558917555928200^^'' 

For N = 1 [ the spin 1/2 Ising model], we have 

X4(l,/3) = -2 - 48/3 - 636/3^ - 6464/3^ - 55892/3* - 2174432/5/3^ - 47009464/15/3** - 2239468288/105/3' - 14570710772/105/3' 
-823130010272/945/3® - 25080975789304/4725/3" - 1640401398782848/51975/3" - 28654566671774104/155925/3^^ 
-2130434175575247424/2027025/3^'' - 83969257269976828688/14189175/3" - 6995762565293277161216/212837625/3^^ 
-38389375874347206695732/212837625/3^® - 272537955948789968719904/278326125/3".. 

For N = 2 [the XY model] we have 

X4(2,/3) = -3/2 - 18/3 - 963/8/3^ - 1233/2/3^ - 171687/64/3* - 167661/16/3^ - 38749413/1024/3" - 32973957/256/3' 
-2142639141/5120/3* - 13411622379/10240/3® - 3907085119879/983040/3" - 240713424017/20480/3^^ 
-1247905418479081/36700160/3^^^ - 166057186013983/1720320/3" 
-407002859073704999/1509949440/3^* - 1960425200264079271/2642411520/3^^ 
-3835682132124206551811/1902536294400/3" - 245360122597207497559/45298483200/3^'... 

For N = 3 [the Heisenberg classical model] , we have 

X4(3,/3) = -6/5 - 48/5/3 - 1076/25/3^ - 11072/75,8'' - 677044/1575/3* - 981856/875/3^ - 958584296/354375/3" 
-261075968/42525/3' - 362572843588/27286875/3* - 9268612328224/334884375/3® 
-3647348945492264/65302453125/3" - 1655707479102099328/15084866671875/3" 
-3670221101428789064/17405615390625/3^^ - 53935317813474946624/135763800046875/3" 
-99068754350666844524336/134632435046484375/3" - 4880947680478330092600064/3635075746255078125/3^'^ 
-192879202499123356385626829692/79771737251567689453125/3" 
-38137398242459901685390609504/8863526361285298828125/3"... 



19 



APPENDIX B: THE SECOND FIELD DERIVATIVE OF THE SUSCEPTIBILITY ON THE BCC 

LATTICE 

The HT expansion coefficients of the second field derivative of the susceptibiHty on the bcc lattice are 

-64 

-2304 - 1200 JV 

''^W = ^3 (2 + ^) 

, , , -32256 - 17408 N 

. , , -3086848 - 4038528 N - 1679424 AT^ - 215600 iV^ 

dAN\ = T 

^ ' N^(2 + Nf (4 + AT) 

-33383424 - 44140288 N - 18541952 - 2396160 
~ JV6 (2 + Nf (4 + N) 

-4025769984 - 8048769024 JV - 6183804416 - 2256738944 - 386061056 JV" - 24592512 

N'' (2 + Nf (4 + AT) (6 + N) 



deiN) 



For the coefficients which follow it is typographically more convenient to set dr{N) = Pr{N)/Qr{N) and to tabulate 
separately the numerator polynomial Pr{N) and the denominator polynomial Qr{N), 

Pt{N) = -38338166784 - 76883050496 N - 59262616576 - 21695712768 - 3720514560 N'^ - 237404160 

QjiN) = JV^ (2 + Nf (4 + JV) (6 + JV) 
PsiN) = -22398548705280 - 64579066675200 Af - 78988001214464 JV^ - 53316161983488 AT^- 



QsiN) = AT^ (2 + AT)* (4 + Nf (6 + AT) (8 + AT) 



P9(JV) = -197760671023104 - 570072979439616 JV - 697267355156480 JV"' - 470715808290816 JV" 



Q9(JV) = (2 + AT)* (4 + JV)^ (6 + JV) (8 + JV) 



Pio(JV) = -135893648670720000 - 506812733453762560 JV - 837961375615287296 AJ"" 



Qio(N) = AJ" (2 + AT)^ (4 + Nf (6 + AT) (8 + AT) (10 + AJ) 



Pii(JV) = -1141044912712581120 - 4249040663391240192 JV - 7015555079356284928 JV"' 
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-6764393828229578752 JV^ - 4222404584768831488 AT* - 1787577880615559168 
-522623903804669952 At" - 105313027917844480 AT^ - 14308686372335616 
-1246403878745600 JV^ - 62559888474112 JV^° - 1369193398272 AT^^ 
Qii(Ar) = AT^^ (2 + AT)-' (4 + AT)'' (6 + N) (8 + N) (10 + N) 
Pi2{N) = -1353224519826548981760 - 6043648174007450075136 AT - 12246178519860941684736 
-14898334498567895384064 Af^ - 12137516326724153901056 AT* - 6993670039948027871232 
-2934699277523969310720 at'' - 909774043355138494464 AT - 209115782694899085312 AT* 
-35427294974940029952 N'-' - 4350609416485181440 Af^" - 375257673683951872 N^^ 
-21477632431608320 AT^^ - 729782591023616 AT^^ - 11103491816320 AT^* 
Qi2iN) = at" (2 + Nf (4 + Nf (6 + Nf (8 + AT) (10 + N) (12 + N) 
Pi3{N) = -10962814317012556185600 - 48855640313305403228160 AT - 98790788652478570168320 
-119949033747113622634496 AT^ - 97539319548016125280256 Af* - 56104431660229407899648 A?*^ 
-23504510497602793209856 Af" - 7275637584279889346560 AT'' - 1670046301220816199680 AT* 
-282577439283542464512 A?" - 34662501972894392320 Ar^° - 2986763073027756032 Af" 
-170792230703923200 AT^^ - 5798723600550656 AT^^ - 88165648564224 AT" 
Qi3{N) = at" (2 + Nf (4 + Nf (6 + Nf (8 + AT) (10 + AT) (12 + AT) 
Pi4(Ar) = -58832199553811357682892800 - 319656261223016276133150720 AT 
-804506884372627522672656384 AT^ - 1245051701802465742955741184 AT^ 
-1326506792510663276202295296 N* - 1032281846527637618727845888 
-607550530706546656156057600 AT^ - 276285749123163236534124544 AT'' 
-98336895868270792502870016 - 27578892092469750017228800 
-6104887252432902432088064 Af^° - 1063450425127326143963136 Af" 
-144612518832699900119040 AT^^ - 15134717223692504156160 AT^^* 
-1191613756119982431744 at" - 68071027536755782912 AT "^'^ 
-2655538734547345536 AT^** - 63104921069570176 AT" - 687220401024000 AT^* 
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Qi4(JV) = JV^^ (2 + (4 + AT)* (6 + JV)^ (8 + JV) (10 + N) (12 + JV) (14 + N) 
(AT) = -463996271577436867233054720 - 2515152883396655909638766592 JV - 6315643027876160842940547072 
-9752398559014235653061738496 - 10368212405864877478064947200 AT* 
-8051953080310134001281007616 - 4729721489677266022278627328 
-2146872892026044965004247040 JV'' - 762790924816930341473615872 JV^ 
-213577145756373531731427328 - 47205533784466451589431296 
-8211418569550233134678016 AT^^ - 1115168266958340559364096 AT^^ 
-116570698649671407411200 Af^^ - 9168049245864669396992 Af" 
-523207897673566349312 AT^^ - 20392833008810153984 AT^** - 484218500284702720 N" - 5269451094097920 AT^® 
Qi5(N) = AT^'^ (2 + AT)'' (4 + AT)" (6 + Nf (8 + N) (10 + N) (12 + N) (14 + N) 
Pie{N) = -3706032155932253013818065551360 - 23511532957902177899390315790336 AT 
-69984745734719724031602453381120 - 129942518358149653545095685734400 
-168777163612514463301490677121024 N* - 163020601132834974169821018587136 
-121523005151702351276069798019072 Af" - 71619975067826655255089882595328 
-33913954817349955356957471670272 AT^ - 13042529099485925377885106077696 n" 
-4101432698345456358253122813952 AT^" - 1058465798499513773553711382528 Af" 
-224369655134353441277932077056 AT^^ - 38996967102960612063728975872 AT^^ 
-5532347732173533521822597120 At" - 635702582117166696141393920 AT^*^ 
-58478148624089490157039616 N^" - 4233637928944288589514240 n" 
-235251240860669046267648 Af^** - 9658642370970465608960 Af" 
-275363031292904868480 - 4856360702431679584 N^^ - 39817256443223280 N^^ 
Qie(Ar) = AT^^ (2 + Nf (4 + Nf (6 + Nf (8 + Nf (10 + N) (12 + AT) (14 + AT) (16 + N) 
Pij{N) = -28625170809317622639273757900800 - 181166520457100717483475767132160 AT 
-537991814958099333262399103827968 AT^ - 996597299051903297408062545985536 N^ 
-1291524895757272450924715899355136 AT* - 1244744825531489741899730562056192 N^ 
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-925926421634294684791564008423424 JV" - 544586262055545477273422842036224 JV'' 
-257372711584126360801923714514944 JV^ - 98795363631395088253841625841664 JV^ 
-31012859792452238511372108300288 JV^° - 7990178234700265757218112274432 AT^^ 
-1691062951010009264622059913216 JV^^ - 293484614468719284469043789824 N^^ 
-41578182122239994876120399872 AT" - 4771486216282518354479939584 JV^^ 
-438407834252226343769522176 iV^" - 31704698511771193169897472 
-1759972842290831947573248 Af^** - 72192427708156325499904 Af^" 
-2056451245211192506368 Ar^° - 36240554109317191808 Af^^ - 296932880108021760 AT^^ 
Qn{N) = AT^* (2 + AT)* (4 + Nf (6 + Nf (8 + AT)^ (10 + AT) (12 + AT) (14 + AT) (16 + AT) 

In particular for N = 0,we have (in terms of the variable f3 = f3/N) 

Xi0) = -3 - 96/3 - 1728, - 24192/3^ - 289392^" - 3129696/3^ - 31451328^® - 299516928^'' - 2734197840/3* 
-24140706912/3® - 207357252000;9^° - 1741096363392/3" - 14339282987664/3^^ - 116166160551840/3^^ 
-927691099407360;9" - 7316490196952064;3^^ - 57068641022992368/3^** - 440794771444139040;9^'' 

For A'' = 1 [ the spin 1/2 Ising model], we have 

X4(l,/3) = -2 - 64/3 - 1168/3^ - 49664/3/3^ - 601360/3/3* - 32820608/15/3^ - 996463616/45/3^ - 66712488448/315/3'' 
-122056132496/63/3* - 48489867797888/2835/3® - 2078558044733696/14175/3^° 
-191285725186144768/155925/3" - 188087379936809600/18711/3^^ - 492034524872707515136/6081075/3" 
-27296494302637993572352/42567525/3" - 3201197677867739316248576/638512875/3^^ 
-4945781553886665074906384/127702575/3^" - 3212941768987291424807915648/10854718875/3" 

For A/' = 2 [ the XY model], we have 

X4(2,|8) = -3/2 - 24/3 - 441/2/3^ - 1572/3^ - 153175/16/3* - 104981/2/3^ - 68238647/256/3° 
-40884131/32/3'' - 1404217891/240/3* - 66125311269/2560/3® - 16313466298147/147456/3^° 
-85154694896333/184320/3" - 51968444431571323/27525120/3^^ - 156353752523792639/20643840/3^ 

-236982809408746803649/7927234560/3" - 400259785750849937/3440640/3^^ 
-637876955227851294690449/1426902220800/3"^" - 48476858619011264835409/28538044416/3"... 

For N = 3 [the Hciscnbcrg classical model], wc have 

X4(3,/3) = -6/5 - 64/5/3 - 1968/25/3^ - 5632/15/3'' - 36138448/23625/3* - 132459392/23625/3^ - 1348186624/70875/3" 
-194131778048/3189375/3' - 228357648983536/1227909375/3* - 2016513048715136/3683728125/3® 
-3366834959446902016/2154980953125/3" - 1872724144398680576/430996190625/3" 
-2674610910995288182912/226273000078125,8" - 367131807235133274368/11636897146875/3" 
-100350171877191930199645952/1211691915418359375/3" - 778819344453674012213276672/3635075746255078125/3" 
-130937251648404408278387043967856/239315211754703068359375/3" 
-22035313880683966424477522278528/15954347450313537890625/3" 



23 



